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We examine the statics and dynamics of vortices in the presence of a periodic quasi-one dimen¬ 
sional substrate, focusing on the limit where the vortex lattice constant is smaller than the substrate 
lattice period. As a function of the substrate strength and filling factor, within the pinned state 
we observe a series of order-disorder transitions associated with buckling phenomena in which the 
number of vortex rows that fit between neighboring substrate maxima increases. These transi¬ 
tions coincide with steps in the depinning threshold, jumps in the density of topological defects, and 
changes in the structure factor. At the buckling transition the vortices are disordered, while between 
the buckling transitions the vortices form a variety of crystalline and partially ordered states. In the 
weak substrate limit, the buckling transitions are absent and the vortices form an ordered hexagonal 
lattice that undergoes changes in its orientation with respect to the substrate as a function of vortex 
density. At intermediate substrate strengths, certain ordered states appear that are correlated with 
peaks in the depinning force. Under an applied drive the system exhibits a rich variety of distinct 
dynamical phases, including plastic flow, a density-modulated moving crystal, and moving floating 
solid phases. We also find a dynamic smectic-to-smectic transition in which the smectic ordering 
changes from being aligned with the substrate to being aligned with the external drive. The different 
dynamical phases can be characterized using velocity histograms and the structure factor. We dis¬ 
cuss how these results are related to recent experiments on vortex ordering on quasi-one-dimensional 
periodic modulated substrates. Our results should also be relevant for other types of systems such 
as ions, colloids, or Wigner crystals interacting with periodic quasi-one-dimensional substrates. 

PACS numbers: 74.25.Wx,74.25.Uv,74.25.Ha 


I. INTRODUCTION 

Commensurate-incommensurate transitions are rele¬ 
vant to a number of condensed matter systems that can 
be effectively described as a lattice of particles inter¬ 
acting with an underlying periodic substrate. A com¬ 
mensurate state occurs when certain length scales of 
the particle lattice match the periodicity of the under¬ 
lying substrate, such as when the number of particles 
is equal to the number of substrate minimaii^. Typi¬ 
cally when commensurate conditions are met, the sys¬ 
tem forms an ordered state free of topological defects, 
while at incommensurate fillings there are several pos¬ 
sibilities depending on the strength of the substrate. If 
the substrate potential is weak, the particles maintain 
their intrinsic lattice structure which floats on top of the 
substrate, while for strong substrates a portion of the 
particles lock into a configuration that is commensurate 
with the substrate while the remaining particles form 
excitations such as kinks, vacancies, or domain walls. 
At intermediate substrate strengths, the lattice order¬ 
ing can be preserved but there can be periodic distor¬ 
tions or rotations of the particle lattice with respect to 
the substrate lattice^""—. These different cases are as¬ 
sociated with differing dynamical responses of the par¬ 
ticles under the application of an external drive^iir— . 
When kinks or domain walls are present, multi-step de¬ 
pinning process can occur when the kinks become mo¬ 
bile at a lower drive than the commensurate portions 
of the sampleii^. Examples of systems that exhibit 


commensurate-incommensurate phases include atoms ad¬ 
sorbed on atomic surfacesi^^i^, vortices in type-II super¬ 
conductors interacting with artificial pinning arraysi^Ti^l, 
vortex states in Josephson-junction array a^^i^^ , super- 
fluid vortices in Bose-Einstein condensates in the pres¬ 
ence of co-rotating optical trap arrays2i“i2&, cold atoms 
and ions on ordered substrates2i“i^, and colloidal par¬ 
ticles on periodic)^-— and quasi-periodic optical 
substrates^i^. 

In the superconducting vortex system, commensura- 
bility occurs when the number of vortices is an inte¬ 
ger multiple of the number of pinning sites, and various 
types of commensurate vortex crystalline states can oc¬ 
cur with different symmetrie a^^i^^i^^i^^ . At fillings where 
there are more vortices than pinning sites, it is pos¬ 
sible to have multi-quantized vortices occupy the pin¬ 
ning sites, and a composite vortex lattice can form that 
is comprised of individual or multiple flux-quanta vor¬ 
tices localized on pinning sites coexisting with vortices 
located in the interstitial regions between the pinning 
sitesi^ii^i22. Ordered commensurate vortex states have 
been directly imaged with Lorentz microscopyi^ and 
other imaging techniques^i^, and the existence of corn- 
mensuration can also be deduced from changes in the 
depinning force needed to move the vortices when peaks 
or steps appear in the critical current as a function of 
vortex densit 3 *i^ii^“— It is also possible for ordered 
vortex structures such as checkerboard states to form at 
rational fractional commensuration ratios of n/m with 
integer m and n, where n is the number of vortices and 
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m is the number of pinning sites^“— . Experiments^ and 
simulations^!^ of colloidal assemblies on optical trap ar¬ 
rays examined the depinning transitions and subsequent 
sliding of the colloids and show that the depinning thresh¬ 
old is maximum for one-to-one matching of colloids and 
traps, while it drops at incommensurate fillings due to 
the presence of highly mobile kinks, anti-kinks, and do¬ 
main walls. 

Commensurate-incommensurate transitions can also 
occur for particles interacting with a periodic quasi-one- 
dimensional (qlD) or washboard potential, where the 
particles can slide freely along one direction of the sub¬ 
strate but not the other. An example of this type of 
system is shown in Fig. [T] for a two-dimensional system 
of vortices interacting with a quasi-one dimensional si¬ 
nusoidal substrate. The potential maxima are indicated 
by the darker shadings, and the vortices are attracted 
to the light colored regions. This type of system has 
been studied previously for colloids interacting with qlD 
periodic substrate arrays, where it was shown that var¬ 
ious melting and structural transitions between hexago¬ 
nal, smectic, and disordered colloidal arrangements can 
occur—-—. In general, the colloidal studies focused on 
the case where the particle lattice constant a is larger 
than the substrate lattice constant w. Martinoli et al 
investigated vortex pinning in samples with a ID peri¬ 
odic thickness modulation^-— and observed broad corn- 
mensuration peaks in the depinning threshold that were 
argued to be correlated with the formation of ordered 
vortex arrangements that could align with the substrate 
periodicity. Other vortex studies in similar samples 
also revealed peaks in the critical depinning force as¬ 
sociated with commensuration effects^!^, while studies 
of vortices interacting with ID magnetic strips showed 
that commensurate conditions were marked by depinning 
steps rather than peaks^. Under an applied dc drive, de¬ 
pinning transitions occur into a sliding state, and when 
an additional ac drive is added to the dc drive, a series of 
Shapiro steps in the voltage-current curves appears when 
the frequency of the oscillatory motion of the vortex lat¬ 
tice over the periodic substrate locks with the ac driving 
frequency^. Similar commensuration effects and Shapiro 
step phenomena were also studied for vortices interact¬ 
ing with periodic washboard potentials or qlD periodic 
sawtooth substrates^. Vortices interacting with periodic 
qlD planar defects have also been studied in layered su¬ 
perconductors when the field is aligned parallel to the 
layer directions. Here, different vortex lattice structures, 
smectic states, and oscillations in the critical current oc¬ 
cur as a function of applied magnetic field^-— . 

For higher vortex densities in the presence of a qlD 
substrate where the vortex lattice constant a is smaller 
than the substrate lattice constant, a < w, there are sev¬ 
eral possibilities for how the vortices can order. In the 
weak substrate limit, they can form a hexagonal lattice 
containing only small distortions, while in the strong sub¬ 
strate limit the vortices can be strongly confined in each 
potential minimum to form ID rows, so that the over¬ 





FIG. 1: The real space images (left column), with the sub¬ 
strate minima indicated by lighter regions and the vortex posi¬ 
tions marked with circles, and the structure factor S{k) (right 
column), for a system with a periodic quasi-one-dimensional 
substrate with Fp = 1.5. (a,b) At w/a = 1.85, each substrate 
minimum contains a single row of vortices (ri) and the struc¬ 
ture factor shows smectic order. (c,d) wja — 2.054, at the 
onset of a buckling transition. (e,f) At wja — 2.651 there is 
an ordered zig-zag r 2 vortex lattice. (g,h) At wja — 3.0 there 
is a mixture of and ra lattices. 


all two-dimensional (2D) vortex structure is anisotropic. 
Between these limits, the vortices can exhibit buckling 
transitions by forming zig-zag patterns within individual 
potential troughs, so that for increasing vortex density 
there could be a series of transitions at which increas¬ 
ing numbers of rows of vortices appear in the potential 
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troughs. Transitions from ID rows of particles to zig¬ 
zag states or multiple rows have been studied for par¬ 
ticles in single qlD trough potentials in the context of 
vortices^-—, Wigner crystals^”—, colloid s qlD 
dusty plasma s ions in qlD traps^“— , and other 
systems^i^ where numerous structural transitions, dif¬ 
fusion behavior and dynamics can occur. In the case of a 
periodic array of channels such as shown in Fig. [U much 
less is known about what buckling transitions would oc¬ 
cur and what the dynamics would be under an applied 
driving force. Recently Guillamon et al. studied vortex 
lattices in samples with a periodic qlD array of grooves. 
As a function of the commensuration ratio p = wja, they 
found that for p < 6, the vortex lattice remains triangular 
but undergoes a series of transitions that are marked by 
rotations of the angle 0 made by the vortex lattice with 
respect to the substrate symmetry direction^S. They also 
observed that at much higher fields, the system transi¬ 
tions into a disordered state with large vortex density 
fluctuations. Open questions include what happens to 
these reorientation transitions as the substrate strength 
is increased, and what the vortex dynamics are when a 
driving force is applied. Dynamical phases and structural 
transitions between different kinds of nonequilibrium vor¬ 
tex flow states have been extensively studied for driven 
vortex systems interacting with randontSl"^— and 2D pe¬ 
riodic pinning arrays^”—; however, there is very little 
work examining the dynamic vortex phases for vortices 
moving over qlD periodic substrates. It is not known 
whether the vortices would undergo dynamical structural 
transitions or exhibit the same types of dynamic phases 
found for vortices driven over random disorder, such as a 
disordered plastic flow state that transitions to a moving 
smectic or anisotropic crystal as a function of increasing 
drive. 

In this work we consider ordering and dynamics of vor¬ 
tices interacting with a periodic qlD sinusoidal poten¬ 
tial for fillings 0<u;/a<5.5. In the strong substrate 
regime, the system undergoes a series of structural tran¬ 
sitions that are related to the number of rows of vor¬ 
tices that fit within each substrate trough. These tran¬ 
sitions include transformations from ID vortex rows to 
zig-zag patterns that gradually increase The vortex 
structure contains numerous dislocations at the buckling 
transitions and is ordered between the buckling transi¬ 
tions. For strong substrates, the onset of the buckling 
and ordered phases produces a series of steps in the crit¬ 
ical depinning threshold as a function of vortex density, 
while for weaker substrate strengths, some of the states 
in which the vortices order produce peaks in the critical 
depinning force. For the weakest substrates, the vortices 
form a triangular lattice that undergoes rotations with re¬ 
spect to the underlying substrate symmetry direction as a 
function of applied magnetic field, similar to the behavior 
observed by Guillamon et al^. Under an applied drive 
we observe plastic flow states, moving density-modulated 
crystals, and dynamic floating solids. For certain fillings 
we also find smectic-to-smectic transitions where the two 


smectic states have different orientations. These different 
flowing phases produce distinct features in the velocity 
histograms and the structure factor. The commensura- 
bility ratio w/a also strongly affects the driving force at 
which the transition to a moving floating solid occurs. 
Our results should be general to other types of systems 
that can be represented as a collection of repulsive par¬ 
ticles interacting with a periodic qlD substrate, such as 
colloids on optical line traps, ions in coupled traps, and 
Wigner crystals on corrugated substrates. 

II. SIMULATION 

We model a two-dimensional system of vortices in¬ 
teracting with a periodic qlD potential with period re, 
where there are periodic boundary conditions in the x 
and ^-directions. The vortices are modeled as point par¬ 
ticles and the dynamics of an individual vortex i obeys 
the following equation of motion: 

V^=Kv+Ft + K + Ffr. ( 1 ) 

Here 77 is the damping constant which we set 
equal to unity. The vortex-vortex forces = 

FoKiiRij/XjKij, where Fq = (/)^/ 27 r/UoA^, (/)o is the 
elementary flux quantum, p is the permittivity, Ki is the 
modified Bessel function, is the location of vortex i, 
Rij = |Ri — Rjl, Rij = (Rf — Tij)/Rij, and A is the 
penetration depth. The vortices have repulsive interac¬ 
tions and form a triangular lattice in the absence of a 
substrate. The vortex interaction with the substrate is 
given by F| = — VU(x^)x where the substrate has the 
sinusoidal form 

V{x) = Vb sin(27rx/u;). (2) 

We define the pinning strength of the substrate to be 
Fp = 27 tVoIw. The dc driving force FJ^ arises from the 
Lorentz force induced by a current applied along the easy 
direction (y-axis) of the substrate which produces a per¬ 
pendicular force on the vortices and causes them to move 
in the x-direction in our system. We measure the vortex 
velocity {Vx) along the driving direction as we increase 
the external drive in increments of and average the 
vortex velocities over a fixed time in order to avoid any 
transient effects. The thermal forces Ft are modeled 
as random Langevin kicks with the properties (Ft) = 0 
and (F^(t)F^(t')) = 2r]kBT6ij6{t — t'), where ks is the 
Boltzmann constant. The initial vortex positions are ob¬ 
tained by annealing from a high temperature state and 
cooling down to T = 0. The dc drive is applied only 
after the annealing procedure is completed. We consider 
a range of vortex densities, which we report in terms of 
the ratio w/a oi the periodicity of the substrate to the 
vortex lattice constant that would appear in the absence 
of a substrate. We denote a state containing n rows of 
vortices in each potential minimum as r^. 
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III. PINNED PHASES 

In Fig. [TJa,c,e,g), we plot the real space locations of 
the vortices on the potential substrate after annealing 
for a system with Fp = 1.5 at fillings of w/a = 1.85, 
2.054, 2.651, and 3.0, while in Fig. [Hb,d,f,h) we show 
the corresponding structure factors S{k). At re/a = 1.58 
in Fig.llja), the vortices form single ID rows in each po¬ 
tential minimum, corresponding to an ri state, and the 
overall vortex structure is highly anisotropic with lattice 
constants ax = 4.5 in the x—direction and ay = 1.31 
in the ^—direction. Additionally, each potential trough 
captures a slightly different number of vortices, introduc¬ 
ing disorder in the alignment of rows in adjacent minima, 
and leaving the system with periodic ordering only along 
the x-direction. The corresponding structure factor in 
Fig. d^b) exhibits a series of spots at ky = 0 and finite 
kx^ indicative of the ID ordering associated with a smec¬ 
tic phase. As the magnetic field increases, the vortex 
ordering must become increasingly anisotropic in order 
to maintain single rows of particles in each minimum. 
This is energetically unfavorable, so instead a transition 
occurs to a zig-zag or buckled state in which there are 
two partial rows of vortices in every substrate minimum. 
Figure [TJc) illustrates the real space vortex positions for 
wja = 2.054 at the beginning of the zig-zag transition, 
where some of the troughs contain a buckled vortex pat¬ 
tern. In the corresponding S{k) plot in Fig. [Hd), the 
smectic ordering develops additional features at large k 
associated with the shorter range structure that arises on 
the length scale associated with the zig-zag pattern. As 
the magnetic field is further increased, the zig-zag pattern 
appears in all the substrate minima and the system forms 
an ordered anisotropic 2D r 2 lattice as shown in Fig.[T]^e) 
for w/a = 2.651, where there are two rows of vortices in 
each potential minimum that form a zig-zag structure 
which is aligned with zig-zag structures in neighboring 
minima. The corresponding S{k) in Fig. [T]^f) has a se¬ 
ries of peaks at small and large k indicating the presence 
of a more ordered vortex structure. For higher fields, 
the zig-zag lattice becomes increasingly anisotropic un¬ 
til another buckling transition occurs to produce rs with 
three vortex rows per substrate minimum. Figure [TJg) 
shows the transition point at w/a = 3.0 where certain 
potential troughs contain three vortex rows while others 
contain two vortex rows or mixtures of two and three 
vortex rows. In Fig. [Hh), S{k) for this case shows that 
the system is considerably more disordered than at the 
commensurate case illustrated in Fig. [He,f). 

In Fig. [2] we show the continuation of the evolution of 
the vortex lattice from Fig. [T]in both real space and k- 
space. At w/a = 3.3 in Fig. [2fa), there is an ordered rs 
structure with three vortex rows in each potential mini¬ 
mum, producing the ordered S{k) shown in Fig.[2fb). As 
the vortex density is further increased, the row structure 
disorders as shown in Fig. He) for w/a = 3.644, corre¬ 
sponding to a ring like structure in S{k) as indicated in 
Fig. [2fd). There are still peaks along the ky = 0.0 line 



FIG. 2: The continuation of the real space images (left col¬ 
umn) and S{k) (right column) from the system in Fig.[T]with 
Fp = 1.5. (a,b) At w/a = 3.3, there is an ordered structure 
with three vortex rows per potential minimum (rs). (c,d) 
At w/a = 3.644, there is a partially disordered state with 
roughly three vortex rows per potential minimum. (e,f) At 
w/a — 4.1455 there is a partially ordered state with four vor¬ 
tex rows per minimum (g,h) The disordered state at 

w/a = 5.15 showing ring structures in S{k). 


due to the anisotropy induced by the substrate. For this 
value of Fp, further increasing the vortex density does 
not produce a more ordered configuration; however, cer¬ 
tain partially ordered structures can occur as illustrated 
in Fig.[2l^e) for w/a = 4.1455, where there are four vortex 
rows per trough (r^) with mixed peaks and smearing in 
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FIG. 3: The fraction of six-fold coordinated vortices Pq vs 
w/a for (a) Fp = 0.5, (b) Fp = 1.0, (c) Fp = 1.5 and (d) 
Fp = 2.0. In panel (c), the labels a, c, e, g indicate the values 
of w/a at which the images in Fig. [T] were obtained, while the 
labels A, C, E, and G indicate the values of w/a at which the 
images in Fig. [2] were obtained. The dips in Pq coincide with 
transitions in the number of vortex rows contained within 
each potential minimum. 


the corresponding structure factor shown in Fig.[2jf). At 
higher fields, the vortex structures become disordered as 
shown in Fig.[2fg) at re/a = 5.15, where S{k) in Fig.[2jh) 
has pronounced ring structures. There are still two peaks 
8it ky = 0 and finite kx due to the smectic ordering im¬ 
posed by the qlD substrate. 

We can also characterize the system using the fraction 
of six-fold coordinated vortices Pq = N~^ 
where Zi is the coordination number of vortex i obtained 
from a Voronoi construction. In general, we find that Pq 
drops at the buckling transitions due to the formation 
of dislocations that are associated with the splitting of 
a single row of vortices into two rows, creating a kink 
at the intersection of the two rows. In Fig. [3l^a) we plot 
Pq versus w/a for a system with Fp = 0.5. Over the 
range 1.0 < w/a < 1.7, each pinning trough contains 
an ri state, while the dip in Pq at w/a = 1.77 corre¬ 
sponds to the middle of the buckling transition when 
there is roughly a 50:50 mixture of ri and r 2 . For 
1.85 < w/a < 2.35, the system forms an ordered r 2 state 
similar to that shown in Fig. [D^e), but less anisotropic 
since the weaker substrate compresses the zig-zag struc¬ 
ture less and permits it to be wider. Near w/a = 2.4, 
there is another buckling transition from r 2 to rs and the 
system forms a disordered state similar to that shown in 
Fig. mg). As re/a is further increased, there is a partially 
ordered state near re/a = 3 that is similar to the state in 
Fig. however, due to the weaker substrate strength 
a fully ordered rs state does not form. For re/a > 3.2 the 


system adopts a polycrystalline configuration that be¬ 
comes more ordered at high vortex densities. In Fig. mb), 
we show that a similar set of features associated with 
buckling transitions occurs for a stronger substrate with 
Fp = 1.0; however, in this case the transition from r 2 to 
rs is sharper and a fully ordered three row state appears 
near w/a = 3.0. 

In Fig. m^) we plot Pq versus w/a for samples with 
Fp = 1.5, the same pinning strength at which the im¬ 
ages in Figs. [Hand [2] were obtained. Here the dips in Pq 
associated with the ri to r 2 , r 2 to rs, and rs to r 4 tran¬ 
sitions are sharper. We also observe the development of 
a small dip near w/a = 4.4 corresponding to a partial 
transition from r 4 to vq . The values of re/a at which row 
transitions occur shift upward with increasing Fp. For 
example, the ri to r 2 transition occurs at w/a = 1.768 
for Fp = 0.5 but at w/a = 2.05 for Fp = 1.5, since the 
higher Fp stabilizes the ri state up to higher lattice con¬ 
stant anisotropies. Figure shows Pq versus w/a for 
samples with Fp = 2.0. Here the dip in Pq at the ri to r 2 
transition broadens, while a pronounced jump emerges 
at w/a = 4.7 corresponding to the r 4 to vq transition. 
We expect that for higher values of Fp, additional dips 
in Pq for transitions from to r^+i states for n > 5 will 
appear ai w/a values higher than those we consider here. 

It is difficult to determine if the buckling transitions 
are first or second order in nature. For particles in an 
isolated trough, the transition from a single row to a 
zig-zag pattern is second order, and there have been sev¬ 
eral studies in cold ion systems of quenches through this 
transition in which the density of kinks was calculated for 
different quench rates and compared to predictions from 
nonequilibrium physics on quenches through continuous 
phase transitions^^-— . We expect that the buckling tran¬ 
sitions we observe are second order; however, it may be 
possible that the additional coupling to particles in neigh¬ 
boring potential minima could change the nature of the 
transition, and we have observed a coexistence of chain 
states which is suggestive of phase separation. For vor¬ 
tex systems it could be difficult to change the substrate 
strength as a function of time, but for colloidal systems 
it is possible to create qlD periodic optical substrates of 
adjustable depth and use them to study time dependent 
transitions by counting the number of kinks that form as 
a function of the rate at which the substrate strength is 
changed. 

In Fig. [4] we plot the depinning force Fc versus w/a 
for Fp = 0.1, 0.25, 0.5, 1.0, 1.5, and 2.0 to show that 
the buckling transitions are associated with changes in 
the slope of the depinning force, which decreases with 
increasing w/a in a series of steps. The first drop in Fc 
near w/a = 2.0 corresponds to the ri to r 2 transition. 
In the ri state, the particle-particle interactions roughly 
cancel in the x-direction, so the depinning force is approx¬ 
imately equal to Fp , while close to the buckling transition 
the vortices on the right side of a zig-zag experience an 
additional repulsive force in the driving direction from 
the vortices on the left side of the zig-zag, decreasing the 
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FIG. 4: The depinning force Fc vs w/a for Fp = 0.1 (red 
circles), 0.25 (green squares), 0.5 (blue diamonds), 1.0 (orange 
up triangles), 1.5 (purple left triangles), and 2.0 (pink down 
triangles) showing that for Fp > 0.25 the buckling transitions 
correspond to step features in Fc. The labels a, c, e, g indicate 
the values of w/a at which the images in Fig. [1] were obtained, 
while the labels A, C, and E indicate the values of w/a at 
which the images in Fig. [2] were obtained. Inset: a highlight 
of the main panel illustrates that for weaker pinning, peaks in 
Fc occur, as shown for Tp = 0.1 (red circles) and 0.25 (green 
squares). The peak is associated with the formation of an 
ordered zig-zag lattice similar to that shown in Fig. [TJe). 



FIG. 5: (a) Real space image of the vortex configuration at 
the peak in Fc at Fp = 0.25 and w/a = 1.94 for the system 
shown in the inset of Fig.|4]where an ordered zig-zag structure 
occurs, (b) The corresponding S{k) contains various peaks 
reflecting the ordered nature of the state. 


driving force needed to depin the vortices. To estimate 
the magnitude of this reduction in Fc, we note that the 
average x-direction spacing between vortices in a given 
zig-zag is approximately dx = 2.0. For a zig-zag with a 
30° angle between the two closest neighbors on the other 
side of the chain from each vortex, the vortex-vortex in¬ 
teraction force of Ki{2) produces an additional repulsive 
force of fr = 0.55, giving a value of Fc = Fp — fr that is 
close to the value of Fc = 0.78 observed after the ri to r 2 
step for the Fp = 1.5 system. Similar arguments can be 
made for the magnitudes of the drops in Fc at the higher 


order transitions as well. At w/a = 1.85 in the Fp = 1.5 
sample, Fc already begins to drop below Fp even though 
the pinned configuration shown in Fig.[T]^a) is an ri state. 
This occurs because in this range of w/a, application of 
a finite Fd < Fc induces a slight buckling of the vortices, 
while for w/a < 1.5 the ri rows remain in a ID pinned 
state up to Fd = Fc. The inset of Fig. 0] shows a blowup 
of Fc versus w/a for the weaker pinning cases Fp = 0.25 
and Fp = 0.1. At rc/a = 1.94 there is a peak in Fc for 
the Fp = 0.25 sample coinciding with the formation of 
the long range ordered zig-zag state shown in Fig. [5]^a). 
The corresponding S{k) in Fig. [5Kb) contains a series of 
peaks that reflect the ordered nature of the state, which 
resembles the zig-zag state in Fig. [IKe,f) except that the 
system is more ordered and the zig-zag structure is wider. 
For Fp = 0.1 the zig-zag state transitions into a hexago¬ 
nal lattice and the peak in Fc begins to disappear. Some 
experiments examining vortices in qlD periodic pinning 
structures show that peaks in the critical current occur 
at certain fUlings^i^'— in regimes where the pinning is 
weak, whereas other experiments performed in the strong 
pinning limit reveal more step-like features in the critical 
current. This suggests that the experiments in the strong 
pinning limit are producing buckling transitions^. 

In Fig. [6] we plot representative real space images with 
the matching S{k) for some other substrate strengths to 
highlight other types of ordering we observe. Figure |6Ka) 
shows the real space ordering of the vortices at Fp = 0.5 
and w/a = 4.825, where the vortex lattice is polycrys¬ 
talline and contains regions of triangular ordering with 
different orientations. The corresponding structure fac¬ 
tor in Fig. IlKb) has ring features with some remnant of 
the smectic ordering appearing at smaller values of k. 
In Fig. IHKc), at Fp = 2.0 and w/a = 4.33 an state 
appears, while S{k) in Fig. EKd) has smectic ordering 
features along with additional crystalline ordering sig¬ 
natures due to the ordered arrangement of the particles 
within the troughs. At Fp = 2.0 and w/a = 4.67 in 
Fig. EKe), a new type of ordered structure appears in 
which the vortices can pack more closely by forming al¬ 
ternate regions of rs and states, producing a consider¬ 
able amount of triangular ordering as seen in the plot of 
S{k) in Fig. EKf), where there are sixfold peaks at large 
k and smectic peaks at smaller k. In Figs. |6Kg,h), for 
Fp = 2.0 and w/a = 5.15, a more disordered structure 
appears, with some regions of the sample containing 
or rs states. 


IV. LATTICE ROTATIONS FOR WEAK 
SUBSTRATES 


We have also studied systems with a pinning strength 
of Fp = 0.02. Here, for w/a > 1.77 the vortices form 
a triangular lattice and the features associated with the 
buckling transitions observed in Fig. [3] are lost. In this 
case the vortex lattice can orient at various angles with 
respect to the underlying substrate. In Fig.[7Ka) we show 
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FIG. 6: Real space images (left column), with the substrate 
minima indicated by lighter regions and the vortex positions 
marked with circles, and S{k) (right column). (a,b) At Fp = 
0.5 and wja — 4.825 there is a polycrystalline structure. (c,d) 
At Fp = 2.0 and wja — 4.33 there is a partially ordered 
structure. (e,f) At Fp — 2.0 and wja — 4.67, an ordered 
structure appears. (g,h) At Fp — 2.0 and wja — 5.15 the 
structure is disordered. 



FIG. 7: Real space images of the vortices in a sample with 
Fp — 0.02 showing that the hexagonal vortex lattice can adopt 
various orientations 0 with respect to the substrate, (a) At 
w/a = 2.8, 0 = 27.6°. (b) At w/a = 3.272, 0 = 24°. (c) At 
w/a = 3.53, 0 = 2.9°. (d) At w/a = 3.75, 0 = 13.2°. 


ing but that the vortex lattice was oriented at an angle 
0 ranging from 0 = 0 to 0 = 30° with respect to the 
substrate. In several cases, they found that the system 
locked to specific angles close to 0 = 30°, 0 = 24°, and 
^ = 0°. We find a much larger variation in the orien¬ 
tation of the lattice with respect to the substrate as a 
function of filling than was observed in the experiments, 
which may be due to differences the pinning strength or 
the finite size of our simulations. Our results show that 
for weak pinning, the buckling transitions are lost and 
are replaced with orientational transitions of the vortex 
lattice with respect to the substrate. Another feature we 
observe when the pinning strength is increased is that 
the vortex lattice becomes disordered or polycrystalline. 
Guillamon et al. also observe that at higher fillings the 
vortex lattice becomes disordered; however, in their sys¬ 
tem there are strong random vortex density fluctuations, 
while for our thermally annealed samples the vortex den¬ 
sity at higher fields is generally uniform. 


the real space vortex positions at w/a = 2.8, where the 
vortices form a triangular lattice that is aligned at an an¬ 
gle 0 = 27.6° with respect to the ^-axis. In Fig. [7|^b) at 
w/a = 3.272, 0 = 24°, while in Fig. (TJ^c) at w/a = 3.53, 
0 = 2.9°, and in Fig. at w/a = 3.75, 0 = 13.2°. 
In Ref;^, Guillamon et al. observed experimentally that 
vortices on a qlD substrate retained triangular order- 


V. DYNAMIC PHASES 

In Fig. [8] we plot simultaneously Pq and the aver¬ 
age nearest neighbor spacing dnn versus F^ for a sam¬ 
ple with Fp = 1.5 at varied w/a. Here, we obtain dnn 
by performing a Voronoi tesselation to identify the Zi 
nearest neighbors of particle i, and then take dnn = 
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FIG. 8: The fraction of six-fold coordinated vortices Pq (dark 
blue curves) and the average nearest-neighbor distance dnn 
(light red curves) vs Fa for a system with Fp = 1.5. (a) 

At w/a = 1.767 the system depins from an ri state, (b) 
At w/a = 2.5 the system dynamically orders into a moving 
triangular lattice with Pq = 0.97. (c) w/a = 3.061. (d) 

At w/a = 3.535, the onset of the dynamically ordered phase 
coincides with a drop in dnn near Fd = 4.0. The labels a, c, 
and e correspond to the values of Fd used for the images in 
Fig. [9l 


ZjLi Gi, where Tij is the distance be¬ 
tween particle i and its jth nearest neighbor. For w/a = 
1.767 in Fig. [8fa), Pq = 0.83 in the pinned ri state that 
occurs for 0 < < 1.4. There is a dip in Pq over the 

range 1.4 < Fd < 1.7, corresponding to the plastic flow 
state in which some of the vortices remain immobile while 
other vortices hop in and out of the potential wells. For 
Fd > 1.7, Pq increases and reaches a saturated value of 
Pq = 0.93 when the vortices form a moving triangular 
lattice containing a small fraction of dislocations. The 
value of dnn drops at the depinning transition, and sev¬ 
eral additional drops in dnn occur at higher drives. In 
the ri pinned state, each vortex has two close nearest 
neighbors that are in the same pinning trough, and four 
more distant nearest neighbors that are in adjacent pin¬ 
ning troughs. Once the vortices depin and enter a moving 
state, they adopt a more isotropic structure, causing dnn 
to drop as the distance to the four more distant near¬ 
est neighbors decreases. The additional drops in dnn ^.t 
higher Fd occur whenever the vortex lattice rearranges 
to become still more isotropic. 

In Fig.[8fb) we plot Pq and dnn versus Fd for the same 
sample at re/a = 2.5 where an ordered zig-zag state with 
Pq = 0.89 appears at zero drive, similar to that shown 
in Fig.lIJe). The depinning threshold is Fc = 0.65, much 
lower than the value of Fc for the w/a = 1.767 filling in 
Fig. [8fa), and the depinning transition is marked by a 
drop in Pq to Pq = 0.2. Over the range 0.65 < Fd ^ 2.0, 
the vortices are in a dynamically disordered state, while 



FIG. 9: Real space images (left column), with the substrate 
minima indicated by lighter regions and the vortex positions 
marked with circles, and S{k) (right column) for the dynamic 
system from Fig. [8jd) with Fp = 1.5 and w/a — 3.535 at the 
values of Fd labeled a, c, and e. (a,b) At Fd — 0.5 the sample 
contains pinned vortices coexisting with individual vortices 
that hop from trough to trough. (c,d) At Fd — 3.5, all the 
vortices move together to form a disordered lattice with a 
periodic density modulations. (e,f) At Fd — 7.0 the system 
forms a moving floating triangular lattice. 


the saturation of dnn above Fd ~ 2.0 indicates that a 
partially ordered state has formed. The value of Pq does 
not reach a maximum until Fd = 3.9, where Pq ~ 0.97 
and an ordered state appears. For 1.5 < Fd < 3.9, we 
observe a moving density-modulated solid. The value of 
the drive at which the sample reaches a moving tri¬ 
angular lattice state, is higher for w/a = 2.5 than for 
w/a = 1.767, even though the depinning threshold Fc is 
smaller for the w/a = 2.5 system. At w/a = 2.5, dnn 
is initially small and jumps up at the depinning transi¬ 
tion, unlike the decrease in dnn at depinning found in 
Fig[8l(a). Since the vortices in Fig.[8l^b) form a zig-zag r 2 
structure in the pinned state, each vortex has four close 
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FIG. 10: Vortex velocities for the system in Fig. [9] with Fp = 1.5 and w/a = 3.535. Left column: Histogram P(14) of 
the instantaneous vortex velocities in the driving direction Vx- Center column: Histogram P{Vy) of the instantaneous vortex 
velocities in the transverse direction Vy. Right column: Heightheld map of Vy versus 14. (a,b,c) The plastic flow regime at 
Fd = 0.5. (d,e,f) The moving modulated solid regime at Fa — 3.5. (g,h,i) The moving floating solid regime at Fa — 7.0. 


nearest neighbors in the same pinning trough, and two 
more distant nearest neighbors in an adjacent pinning 
trough. This causes dnn to be smaller in the pinned state 
than it was for the ri structure in Fig. [Hl^a), and when 
the vortex lattice becomes more isotropic in the moving 
state, dnn increases rather than decreasing as the two 
halves of each zig-zag structure move further apart. At 
Fd = 4.1 we observe a drop in dnn that coincides with 
a dip in Pq. This feature is associated with a transi¬ 
tion from a density-modulated lattice to a more uniform 
moving floating lattice. 

In Fig. IHl^c), we show Pq and dnn versus Fd 8it w/a = 
3.061 where there is an rs pinned state, similar to that 
illustrated in Fig. [2fa). Here the depinning threshold 
Fc = 0.3, and the system transitions into a moving tri¬ 
angular lattice at F/^ = 1.6, which is somewhat lower 
than the value of F/^ for w/a = 2.5 in Fig. [Hfb). The 
behavior of dnn in Fig. [Sfc) follows a similar pattern as 
in Fig. [Sl^b), with dnn increasing with increasing Fd. We 
plot the same quantities for w/a = 3.535 in Fig. (H^d), 
where the depinning threshold Fc ~ 0.087 and the sys¬ 
tem dynamically orders for > 4.0. There is a small dip 
in dnn at = 4.15 along with a saturation in Pq which is 
correlated with a structural change to a dynamic floating 
lattice. 

In order to characterize the nature of the dynamic vor¬ 


tex structures in the moving states, in Fig.[9]^a,b) we plot 
the real space images and S{k) for the system in Fig.[8l^d) 
at w/a = 3.535 and Fd = 0.5. Here, individual vortices 
jump from one pinning well to the next while a portion 
of the vortices remain immobile in the substrate minima. 
As shown in the plot of S{k)^ the vortex configuration 
is fairly ordered and takes the form of a distorted non- 
triangular structure which causes Pq to be low for this 
value of Fd. For Fd > Fp, the vortices move together so 
there is no plastic motion, and form a distorted lattice 
containing pronounced density modulations as shown in 
Fig. [9l^c,d) for Fd = 3.5. For Fd > 4.0 we find a tran¬ 
sition from the density modulated lattice to a moving 
homogeneous floating triangular lattice which coincides 
with the drop in dnn in Fig. [8fd) and the maximum in 
Pq. Fig.[9je) shows the floating lattice at F^^ = 7.0, where 
as indicated in Fig.[9ff) S{k) contains sixfold peaks that 
are indicative of triangular ordering. The smectic order¬ 
ing induced by the substrate is substantially weaker or 
almost absent at this drive, as shown by the weakness 
of the spots in S{k) at ky = 0, indicating that the sys¬ 
tem has formed a floating solid. We find similar types of 
transitions in the dynamics at other fillings as well. 

We can also characterize the different dynamic states 
in Figs. [8] and [9] by examining histograms of the vortex 
velocities. In Fig. uni: a) we plot the distribution F(I4) of 
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FIG. 11: Real space images (left column), with the substrate 
minima indicated by lighter regions and the vortex positions 
marked with circles, and S{k) (right column) for a system 
with Fp = 0.5 and wja = 1.767. (a,b) The plastic flow phase 
at Fd = 0.3 where there is individual vortex hopping from 
well to well. The S(k) peaks indicate a smectic phase with 
periodic ordering along the x—direction. (c,d) At Fa = 0.6, 
all the vortices are flowing and form chains that are aligned 
in the x-direction. S{k) shows that a new smectic order has 
appeared with periodic ordering along the ^-direction. (e,f) 
At Fd = 3.6 there is a moving floating triangular crystal. 


Vx in the driving direction at = 0.5 for the system in 
Fig.[9]^a,b) with Fp = 1.5 and wja = 3.535. Figure [TOT b) 
shows the transverse velocities P{Vy), while in Fig. flOl c) 
we plot Vy versus 14 as a heightfield map. At this drive, 
the motion is plastic and occurs by individual vortex hop¬ 
ping, so there is a sharp peak in P( 14) at 14 =0.15 which 
reflects the fact that most of the vortices are slowly mov¬ 
ing within an individual pinning trough. When a single 
vortex jumps into an adjacent pinning trough, it creates 
a pulse of motion through the trapped vortices that trig¬ 
gers the jump of another single vortex to the next pinning 
trough, where the process repeats. This depinning cycle 
creates two peaks in P(I4). The peak at low Vx cor¬ 


responds to the motion of a velocity pulse through the 
dense assembly of vortices at the bottom of the pinning 
trough, while the peak at high Vx is produced by in¬ 
dividual vortices escaping over the potential maximum. 
This peak falls near Vx = 1.1, which is larger than P^^, 
reflecting the fact that after an individual vortex passes 
the crest of the substrate maximum, the substrate con¬ 
tributes an additional force term in the driving direction 
as the vortex moves toward the next substrate minimum. 
In this case the maximum force exerted by the pinning 
site is Fp = 1.5 while the driving force is F^ = 0.5, so 
that the maximum possible instantaneous vortex velocity 
would be Vx = 2.0; however, vortex-vortex interactions 
prevent individual vortices from moving this rapidly. In 
Fig. [TOT b). P{Vy) is centered at I 4 = 0 since there is 
no driving force in the transverse direction; however, we 
observe some asymmetry in P{Vy) as well as peaks at 
finite Vy due to the fact that the vortex lattice segments 
inside the pinning troughs are oriented at an angle with 
respect to the substrate symmetry direction, as shown in 
Fig. [9l^a). This asymmetry also appears in the Vy versus 
Vx plot in Fig. [TOT c). which has two prominent features. 
The first is a wide band of Vy values at low Vx that are 
associated with soliton-like pulses moving through the 
dense regions of the vortex clusters, which push vortices 
in both the positive and negative ^-direction. The second 
is the loop shape at larger Vx values which corresponds 
to motion in which the vortices are accelerated or de¬ 
celerated as they pass over the substrate maxima and 
minima. 

In Fig. flQlf d.e.f) we show instantaneous velocity plots 
for the system in Fig. [9jc,d) with F^ = 3.5 where the 
vortices are moving elastically in the density-modulated 
solid phase. Here P(14) in Fig. [TUTd) has peaks at 
Vx = 2.01 and Vx = 4.75 which are smoothly connected 
by finite P(14) values. The shape of this histogram shows 
the velocity imposed by the driving force of Fd = 3.5 
when the substrate forces alternately act with or against 
the driving force. When the substrate force is against 
the drive the velocity is Vx = Fd — Fp = 2.0, while when 
the substrate and driving forces are in the same direc¬ 
tion, Vx = Fd Fp = 5.0, close to the observed values 
of the peaks in P(I4). In Fig. [TOlf e). P{Vy) has two 
peaks close to Vy = 0.1 and Vy = —0.1, indicating that 
there is an oscillatory motion in the ^-direction. This 
effect can be seen more clearly in the Vy versus Vx plot 
in Fig. [TOlf f) which has two symmetric lobes. Since the 
vortices are in a density-modulated lattice, shearing in 
the ^-direction occurs between adjacent density modula¬ 
tions, with one density modulation moving in the pos¬ 
itive ^-direction while the other moves in the negative 
^-direction. 

The velocity plots for the system in Fig. at 

Fd = 7.0 appear in Fig. [IDi;g,h,i). In Fig. HD^g), P(I4) 
has a two-peak feature similar to that in Fig. ITTil di. but 
with peak values at Vx = 5.52 and Vx = 8.4. Figure 
Mh) shows that P{Vy) has a single peak centered at 
14 = 0, while in Fig. unii), there is a single lobe in the 
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FIG. 12: Vortex velocities for the system in Fig. [TT] with Fp = 0.5 and w/a = 1.767. Left column: P{Vx). Center column: 
P{Vy). Right column: Heightheld map of Vy versus 14. (a,b,c) The plastic flow phase at Fa = 0.3 where there is a large peak 
in P{Vx) at 14 = 0 due to the pinned vortices. (d,e,f) The moving smectic phase from Fig. fTTl c.d') at Fd = 0.6. (g,h,i) The 
moving triangular solid phase at Fd = 3.6. 


Vy versus 14 plot. Here the vortices have formed a float¬ 
ing triangular solid, and their motion is close to one¬ 
dimensional along the driving direction. As Fd is further 
increased, the width of the lobe feature in the Vy direc¬ 
tion gradually decreases. We find similar histograms for 
the other fillings in the strong pinning limit for the plastic 
flow, moving modulated solid, and moving floating solid 
regimes. 


A. Smectic to Smectic Transitions 

For Fp < 1.0, we find that a dynamically induced smec¬ 
tic to smectic transition can occur. In Fig. [TTla.bl we 
show the real space and S{k) images for a system with 
Fp = 0.5 and w/a = 1.767 in the plastic flow regime 
where there is a combination of vortices that are trapped 
in the pinning troughs and a smaller amount of vortices 
that hop by jumping from one trough to the next and 
then triggering a jump of another vortex from one trough 
to the next. Here, S{k) indicates that the overall system 
has smectic ordering due to the chain-like structure of 
the vortices within the pinning troughs. In Fig. [T2ja) we 
plot P{Vx) at Fd = 0.3, where there is a peak at 14 = 0 
due to the pinned vortices along with a small bump at 
Vx = 0.6 due to the vortex hopping. Figure [T2fb) shows 


that PiVy) has a maximum at Vy = 0, while in Fig.[T2fc), 
the Vy versus Vx plot is asymmetric in with a peak at 
14 = = 0.0 and a second peak at higher Vx produced 

by the moving vortices. 

In Fig. HHc) we show the real space vortex configu¬ 
ration at Fd = 0.6, which is higher than the maximum 
pinning force of Fp = 0.5. All the vortices are in mo¬ 
tion, but instead of retaining their alignment along the 
^-direction induced by the substrate, they form a chain¬ 
like structure aligned in the x-direction with a slight tilt 
in the positive ^-direction. This alignment in the drive 
direction is more clearly seen in the corresponding S(k) 
in Fig. [TH dl. where the peaks fall along kx = 0, indicat¬ 
ing a smectic phase with ordering along the ^-direction. 
There are some very weak peaks on the ky = 0 axis due to 
the substrate, but overall the vortex structure is a smec¬ 
tic state rotated 90° from the ^-axis. The peaks do not 
fall exactly at = 0 but are at a slight angle, due to the 
channels in Fig. me) being slightly tilted in the positive 
^-direction. In Fig. [I2fd), P{Vx) for this case shows a 
peak at Vx = 0.19, while there is an absence of weight 
in P(I4) at Vx = 0.0, indicating that the vortices are al¬ 
ways in motion. Figure [T^ej shows that P{Vy) peaks at 
= 0.0 and has an overall asymmetry, which also ap¬ 
pears in the Vy versus Vx plot in Fig. [T2jf). The vortex 
channeling occurs when the vortices form effective pairs 
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FIG. 13: (a) , the drive at which the system transitions 

from a density modulated moving crystal to an ordered mov¬ 
ing floating solid, vs le/a at Fp = 0.5 (red circles), 1.0 (green 
squares), 1.5 (orange diamonds), and 2.0 (black triangles), 
(b) Dynamic phase diagram as a function of Fd and le/a for a 
system with Fp = 1.5. P: pinned phase. Plastic: plastic flow 
regime. ML: moving modulated lattice state. MFS: moving 
flowing solid state. Dashed lines are guides to the eye that 
indicate the transition from ri to r 2 , r 2 to rs, and rs to a 
disordered pinned state. 


aligned along the x-direction. In each pair, one vortex 
is slowed by the backward-sloping side of the potential 
trough, while the other vortex is sped up by the forward- 
sloping side of the potential. The faster vortex pushes 
the slower vortex, giving the pair an increased net mo¬ 
tion along the x direction. This pairing effect is visible 
in the real space image in Fig. HIc). 


As Fd further increases, there is a transition to a 
flowing solid phase as shown in Fig. fTTl e.f) for Fd = 
3.6, where S{k) has sixfold ordering. The correspond¬ 
ing P{Vx) in Fig. [T^g) has two peaks, while P{Vy) 
in Fig. [T^ h) has a symmetrical distribution with three 
peaks indicating that there is an oscillation in the vor¬ 
tex orbits in the ^-direction. The plot of Vy versus Vx in 
Fig.fT^h) contains a single lobe similar to that found for 
the moving floating solid in Fig. ttUi). As Fd is increased 
still further, the width of this lobe in the Vy direction de¬ 
creases. The smectic-to-smectic transition is limited to 
the range 1.5 < w/a < 2, in which two vortices can fit 
between adjacent potential maxima in the dynamically 
moving regime. 


VI. DYNAMICAL PHASE DIAGRAM 

For Fp > 0.25, the drive at which the system tran¬ 
sitions from a density modulated moving crystal to an 
ordered moving floating solid shows considerable varia¬ 
tion with re/a, particularly for the larger values of Fp. In 
Fig. [T3f a) we plot determined from the location of a 
feature in Pq, versus w/a for Fp = 0.5, 1.0, 1.5, and 2.0. 
For Fp = 0.5, has a local maximum near w/a = 2.5, 
and then drops for re/a > 3.0. In the pinned phase for 
w/a > 3.0, the system forms a polycrystalline state, and 
in the moving state the grains realign to form a mov¬ 
ing crystalline state. For Fp = 1.0, 1.5, and 2.0, when 
w/a < 1.75 the system depins from a single chain of vor¬ 
tices and can partially form a moving crystal state. When 
re/a is large enough that a pinned zig-zag state forms, the 
moving density-modulated state can persist up to much 
higher drives. For Fp = 1.0 the system forms a pinned 
poly crystalline state for w/a > 4.5 which coincides with 
the drop in for w/a > 4.5. For Fp = 1.5 and 2.0, 
there are local peaks in that correlate with the mov¬ 
ing buckled phases which occur when groups of vortices 
can fit between adjacent pinning maxima as the vortices 
move. This effect is most pronounced for Fp = 2.0. 

In Fig. fldfb) we plot a dynamic phase diagram as a 
function of Fd and w/a for a system with Fp = 1.5. 
Above depinning in the regime where Fd < Fp, the sys¬ 
tem is in a plastic flow state in which there is a coexis¬ 
tence of moving vortices and immobile vortices. In this 
regime the structure factor generally shows disordered 
features. For Fd > Fp, all the vortices are moving and 
the system is either in a modulated lattice (ML) state 
or a moving flowing solid (MFS) state. We find similar 
dynamical phase diagrams for other values of Fp > 1.0, 
while for the weaker substrates, the size of the plastic 
flow region is reduced and the ML phase is replaced 
with a smectic moving state similar to that shown in 

Fig.[nic,d). 


VII. DISCUSSION 

The dynamic phases we observe have certain similar¬ 
ities to the dynamic states observed for vortices moving 
over random pinning arrays in that there can be pinned, 
plastic, and dynamically ordered phases as a function 
of external drive^i"— . There are some differences, in¬ 
cluding the fact that the moving modulated lattice we 
observe does not form a smectic state that is aligned 
in the drive direction, as found for vortices moving over 
random pinning arrays^”—. Future studies might con¬ 
sider combinations of random disorder with periodic dis¬ 
order, which would introduce a competition in the mov¬ 
ing phase between the smectic ordering imposed by the 
substrate and the smectic ordering induced by the drive. 
Additionally, for random pinning arrays simulations in¬ 
dicate that for increasing vortex density, the drive at 
which the transition to the ordered state occurs decreases 
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due to the increase in the strength of the vortex-vortex 
interactions^^. For the qlD substrate, the location of 
the ordering transition shows strong fluctuations due to 
the ability of the moving lattice to become dynamically 
commensurate with the periodicity of the substrate. 

VIII. SUMMARY 

We examine the statics and dynamics of vortices inter¬ 
acting with a periodic quasi-one-dimensional substrate 
in the limit where the vortex lattice spacing is smaller 
than the spacing of the periodic lattice. For weak sub¬ 
strate strengths, we find that the vortices retain hexago¬ 
nal ordering but exhibit numerous rotations with respect 
to the substrate, similar to recent experimental observa¬ 
tions. For stronger substrates there are a series of buck¬ 
ling transitions where the vortices can form anisotropic 
ID chains, zig-zag patterns, and higher order numbers of 
chains within each substrate minimum. At some fillings 
the overall lattice has long range order and becomes par¬ 
tially distorted at the transitions between these states. 
For higher fillings the buckling transitions are lost and 
the system forms a polycrystalline state. We also find 
that the depinning shows a series of step like features 


when the system transitions from a state with n chains 
to a state with n + 1 chains in each substrate minimum, 
and that for weaker pinning there are some cases where 
there is a peak in the depinning force as a function of fill¬ 
ing. For weak substrates, under an applied drive the vor¬ 
tices depin elastically and retain their triangular order¬ 
ing, while for strong substrates the buckled states tran¬ 
sition to a partially disordered flowing state followed by 
various other transitions into moving modulated crystal 
or homogeneous floating moving crystal states. Our re¬ 
sults should also be applicable to other systems of par¬ 
ticles with repulsive interactions in the presence of a pe¬ 
riodic quasi-one dimensional substrate, such as electron 
crystals, colloids, and ions in optical traps. 
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